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-^The  five-layered  composite  as  treated  In  this  work  consists  of 
three  plates  bonded  together  with  adhesives.  The  thickness  of  th< 
adhesives,  though  small.  Is  Included  In  the  analysis  resulting  In 
the  five  layer  geometry.  A dominant  flaw  is  assumed  to  prevail 
at  the  center  of  the  middle  layer  when  the  composite  Is  being 
stretched  uniformly  and  to  overshadow  the  other  types  of  mechani- 
cal damages  such  as  matrix/fiber  debonding,  fiber  breaking,  etc. 
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In  order  to  gain  a knowledge  of  the  failure  mode  of  the  layered 
composite,  the  load  transmission  characteristics  which  depend  upon 
the  geometrical  and  material  constants  must  be  known  through  a 
stress  analysis.  This  is  accomplished  by  appealing  to  an  approxi- 
mate theory  of  laminates  developed  from  variational  principles. 

The  character  of  tfra  solution*  particularly  that  portion  around  the 
crack  front, Tis  three-dimensional  and  is  sufficiently  accurate  for 
predicting  physically  realizable  failure  modest  In  addition  to 
analyze  the  influence  of  the  adhesive  layers, cf^he  stress  distribu- 
tion obtained  tese  may  be  further  applied  in  conjunction  with  a 
failure  criterion  to  determine  whether  de lamination  will  occur 
prior  to  unstable  crack  propagation  or  not.  The  pertinent  param- 
eter, stress  intensity  factor,  is  obtained  for  a variety  of  mate- 
rial and  geometrical  constants. - 
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SUMMARY 

The  five-layered  composite  as  treated  In  this  work  con- 
sists of  three  plates  bonded  together  with  adhesives.  The 
thickness  of  the  adhesives,  though  small.  Is  included  In  the 
analysis  resulting  In  the  five  layer  geometry.  A dominant 
flaw  Is  assumed  to  prevail  at  the  center  of  the  middle  layer 
when  the  composite  Is  being  stretched  uniformly  and  to  overshadow 
the  other  types  of  mechanical  damages  such  as  matrix/fiber 

: 

. 

debonding,  fiber  breaking,  etc. 

In  order  to  gain  a knowledge  of  the  failure  mode  of  the 
layered  composite,  the  load  transmission  characteristics  which 
depend  upon  the  geometrical  and  material  constants 
must  be  known  through  a stress  analysis.  This  is  accomplished 
by  appealing  to  an  approximate  theory  of  laminates  developed 
from  variational  principles.  The  character  of  the  solution, 
particularly  that  portion  around  the  crack  front,  is  three- 
dimensional  and  is  sufficiently  accurate  for  predicting  phys- 
ically realizable  failure  modes.  In  addition  to  analyze  the 
influence  of  the  adhesive  layers,  the  stress  distribution  ob- 
tained here  may  be  further  applied  in  conjunction  with  a fail- 
ure  criterion  to  determine  whether  delamination  will  occur 
prior  to  unstable  crack  propagation  or  not.  The  pertinent 
parameter,  stress  intensity  factor,  is  obtained  for  a variety 
of  material  arid  geometrical  constants. 
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INTRODUCTION 


Advanced  laminate  composites  are  finding  Increasing  appli- 
cations in  load  carrying  aircraft  components.  As  a result, 
there  has  been  a growing  need  for  reliable  methods  of  deter- 
mining strength  and/or  life  expectancy  of  composites.  The 
fracture  mechanics  approach  addresses  the  physical  processes 
of  material  damage  that  actually  occur  in  the  composite.  Ana- 
lytical modeling  of  the  variety  of  failures  in  composites  such 
. 

as  matrix/fiber  debonding,  fiber  breaking,  matrix  cracking,  ply 
de lamination,  free  edge  effects,  etc.,  if  carried  out  Indis- 
criminately, yields  little  or  no  useful  information.  The  other 
extreme  measure  of  simply  testing  specimen  after  specimen  is 
no  better  as  very  little  understanding  could  be  gained. 

For  those  laminates  such  as  graphite/epoxy  plates  that  ex- 
hibit brittle  fracture  behavior,  there  prevails  a distinct 
point  of  instability  on  the  load-deformation  curve;  it  is  pos- 
sible to  identify  the  physical  damage  such  as  a dominant  crack 
in  the  composite  Just  prior  to  load  instability.  Attention 
should  then  be  focused  in  determining  the  size,  shape  and  lo- 
cation of  this  dominant  crack.  The  other  damages  such  as  fiber 
breaking,  matrix/fiber  debonding,  etc.,  that  occurred  in  the 
early  part  of  load  history  should  be  considered  as  dissipated 
energy  that  Is  no  longer  available  at  the  final  stage  of  load 
instability.  This  is  a very  brief  description  of  the  idea  on 
which  fracture  mechanics  is  based  upon. 


In  a previous  report  [13,  the  authors  developed  a method 
for  solving  the  stresses  in  a multi-layered  composite  plate 
containing  a through-the-thickness  crack  at  the  center.  A 
three- layered  laminate  was  solved  numerically  as  an  illustra- 
tive example  for  the  crack  front  stress  intensity  factors.  The 
initial  assumption  of  a crack  cutting  through  all  the  layers, 
however,  may  not  be  realistic.  This  is  because  the  zone  of 
hydrostatic  tension  or  high  dilatation  occurs  only  at  the 
center  of  the  middle  layer  of  the  laminate.  In  fact,  the 
outer  layers  distort  relatively  more  and  are  less  prone  to 
fail  by  cracking  through  the  thickness  coinciding  with  that 
in  the  middle  layer.  There  is  the  possibility  that  interface 
or  interlayer  delamination  accompanied  by  crack  propagation 
in  the  middle  layer  may  occur  first  before  the  outerlayers  be- 
gin to  separate*.  The  load  transmission  characteristics  in  a 
layered  composite  can  be  very  complex  and  will  not  be  under- 
stood unless  effective  stress  analysis  is  performed  and  used 
Jointly  with  a realistic  fracture  criterion. 

As  a first  prerequisite,  this  report  attempts  to  model  the 
laminate  with  a center. crack  in  the  middle  layer  while  the 
outerlayers  are  uncracked.  This  obviously  makes  the  problem 
considerably  more  difficult.  In  addition,  adhesive  layers 


This  has  been  observed  in  a three-layered  plexiglass  plate 
that  was  pulled  apart  with  a crack  in  the  middle  layer.  The 
middle  layer  can  be  separated  completely  by  crack  propagation 
with  the  outerlayers  still  intact.  This  test  was  performed  at 
the  Lehigh  University  Institute  of  Fracture  and  Solid  Mechanics. 
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are  inserted  between  the  adjoining  layers.  Hence,  a three- 
layered plate  will,  in  fact,  contain  five  layers  if  the  two 
adhesive  layers  are  also  included.  Whether  the  adhesive 
layers  will  or  will  not  have  a significant  influence  on  de- 
lamination  and/or  crack  propagation  can  only  be  answered  by 
appealing  to  a stress  and  fracture  analysis. 

An  approximate  theory  of  laminates  similar  to  that  given 
in  £l]  is  used.  The  governing  equations  are  derived  from 
the  complementary  potential  energy  theorem  in  variational  cal- 
culus. Near  the  crack  front,  the  approximate  approach  does 
preserve  the  qualitative  character  of  the  solution  which  is 

« 

inherently  three-dimensional  [2].  What  is  different  is  the 
load  transmission  characteristics  that  are  reflected  through 
the  stress  intensity  parameter.  A numerical  scheme  is  de- 
veloped for  finding  these  parameters  for  various  geometrical 
and  material  constants.  The  results  are  reported  herein. 

The  stress  components  around  the  crack  are  also  calculated 
and  presented  which  are  needed  in  the  study  of  delamination. 
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The  principle  of  minimum  complementary  energy  as  used  in 
calculus  of  variations  is  a standard  procedure  for  developing 
approximate  theories  of  plates  and  has  been  widely  applied  in 
the  field  of  plates  and  shells.  An  approximate  form  of  the 
solution,  say  in  terms  of  the  stress  components,  is  first  as- 
sumed which  are  then  substituted  into  the  complementary  poten- 
tial energy  known  as  a functional.  Minimization  of  this  func- 
tional with  respect  to  the  admissible  variations  of  the  stress 
components  leads  to  a best  approximation  of  the  exact  solution 
in  a weighted  sense. 

The  procedure  to  be  employed  for  laminated  plates  is  to 
assume  that  each  of  the  stress  components  in  the  pth  layer 
can  be  approximated  as  a product  of  a function  of  the  trans- 
verse variable  z multiplied  by  a function  of  the  in-plane 
variables  x and  y,  i.e., 

Cax»ay*Txy^  " fp(z)CSx(x,y),  Sy(x,y),  SxyCx,y)] 

CTxz»TyZJ  * * fpCz)[ZxCx,y),  ZyCx,y)]  Cl) 

P P 

[a  ] - f Cz)[Z.(x,y)] 

P v P 

in  which  prime  designates  differentiation  with  respect  to  the 
variable  z.  The  function  fp(z)  describes  the  stress  variations 
in  the  thickness  direction  through  the  pth  layer  of  the  lami- 
nate . _ 


This  assumed  form  of  the  stress  field  is  substituted  into 
the  complementary  potential  energy  functional.  Minimization 


of  the  functional  with  respect  to  f(z),  S , S , S , Z , Z 

x y xy  x y 

and  Zz  of  each  layer  is  then  carried  out  whereby  the  conditions 
of  equilibrium  are  enforced.  The  choice  of  the  stress  and 
displacement  expressions  should  be  such  that  the  approximate 
solution  satisfies  the  continuity  conditions  across  the  mate- 
rial interfaces  as  closely  as  possible.  For  example,  the  trac- 
tion continuity  conditions  are  satisfied  by  taking 


(2) 


Vz)  " fp+lCz)’  fp(z>  “ fi+i<>>  C3) 

in  which  z corresponds  to  the  respective  values  at  the  inter- 
face between  layers  p and  p+1.  This  assumption  precludes  an 
exact  evaluation  of  the  state  of  affairs  in  those  regions  of 
the  outer layers  that  are  close  to  the  crack.  Nevertheless, 
the  solution  is  expected  to  be  accurate  at  distances  sufficiently 
far  away  from  the  crack.  Moreover,  the  qualitative  features  of 
the  crack  front  stresses  will  not  be  altered  by  assumptions  (2) 
and  (3).  The  same  argument  leads  to  the  following  displacement 
continuity  conditions: 


—6 


1 » 


Finally,  the  assumed  stress  field  takes  the  form 


(5) 

/ 

The  stress  equilibrium  conditions  in  three-dimensional 

elasticity  when  enforced  on  the  stresses  a , a , , etc.,  in 

x y 

equations  (5)  yield  a set  of  equations  involving  the  variables 
x and  y only: 


(6) 


Note  that  once  S , S and  S are  known,  the  quantities  Z , 

X J AJ  A 

Zy  and  Zz  that  govern  the  transverse  stress  distribution  can 
be  obtained  by  means  of  equations  (6).  This  simplification  is 
a direct  consequence  of  the  product  form  of  solution  assumed 
in  equations  (1). 

The  Va.fu.nt4.onat  Palnclptt  of  Utntmum  Comptementaxy  Pote.nt4.at 

Enixgy  {ok.  Laminate 4.  The  theorem  of  minimum  complementary 

potential  energy  in  variational  calculus  may  be  used  to  develop 
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as  as_ 

X I Aj 

37"  + -ay 


as  as 


3Zx  3Zy 

W + JT 


^■ax,ay*Txy^p  * fp ^Sx,Sy,Sxy^ 

^Txz,Tyz^p  " “ fp^z^Zx,Zy^ 
[a2]p  - fpCz)Zz 


an  approximate  theory  for  laminated  plates.  This  approach  was 
used  In  [1]  and  the  same  procedure  is  applicable  to  the  current 
system.  Without  going  into  the  details  [lj,  the  differential 


equations  governing  the  stresses  S , S , , Z and  the  weighted 

A jf  j 

displacements  u_,  u„  and  u can  be  written  as 
a y z 


3u 


ZlSx  + J2Sy  + hZz  + JT 


3u. 


hSy  * V*  * hzz  * 1y*  m 0 
3u  3u 

2(Il-I2)Sxy  + W"  + 3?  " 0 


I3Zx  + ux  + 3T  9 u 
3uz 

12  + + -JF— • ■ Q 

3 y y 3y 


3u. 


(7) 


WV  + i5zz  + uz  - o 


in  which  the  following  contractions  have  been  made: 
l hl 

I1  " IT  l CfiCz)32dz»  z2  “ “ vixl 
1 ~hl 

2C1+V,)  hl  v,  hl 

I3  « jj*-  / Cf{(z)32dz,  I,,  - - / f^zJfJCzJdz 

Hx  1 1 (8 

Z5  9 h [ Cf1(z)]2dz 
1 "hl 
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In  equations  (7),  the  weighted  displacement  components  ux, 

u and  u are  related  to  the  actual  displacements  0 , fJ  and 
* x y 

!JZ  through 


u » - / ff_f"Cz)dz 


“y  ■ - l V;(z)da 
p 


V£Cz)dz 


where  Hp  denotes  the  thickness  of  the  pth  layer. 

Equations  (6)  may  be  used  to  reduce  the  system  of  equations 
(7)  in  terms  of  the  unknowns  Z%i  Z and  uz.  The  results  are 


z*  - a67’zx  - 51  (al7*uz  + a2uz>  ' 0 


Zy  ‘ a67'Zy  - 5y  <al7*uz  + a2uz>  * 0 


a37X  * V*uz  + a5uz  ' 0 


32  3Z 

53T  ♦ 5y  + VX  + a8uz  * 0 


The  symbols  7*  and  7*  represent  the  harmonic  and  biharmonic 
operators  in  the  variables  x and  y while  a1#  a 2,- — , ag  are 
complicated  functions  of  I^t  I2, , 1^  as  defined  in  equa- 
tions (8),  i.e.. 


I^CI?-2I1|)+2I1II. 

a!  “ 5TTj=rp c i4a3-2i4)+t5ci1+t2 : 


2Ir^ 

:4cr3"-;>i4)+t;ti1+i2: 


a3  - zi  - v5 


a4  - - 21^-12 ) + IXI3 


a5  " Hipip  • a6 


X3 


l4  u3-^x4  ;fi5U2+JL2 . 


aQ  » 


V*2 


.4U3-<'i4;+i5U1+i2, 


Once  Zx>  Zy  and  uz  are  found,  the  remaining  unknowns  ux,  uyt 
etc.,  can  be  obtained  in  a straightforward  fashion  since 


ux  " 


- hZx  ' x 


3u 

“ I3Zy  " W 
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1 , > Uz  8 UZ  8ZI 

s*  ■ ipif  (I1  Ihh  - W1*”  >T 


- t-hh  + WV1  a/> 


C!  , 1 r T 32UZ  , 

sy  iptf  (I1  IF"  ' 1 


2 ^I1I3  - 3^ 


- [I2I3  + X^^-Xj)]  ^S) 


1 »*u* 

S - 5. 

xy  I1-I2  3x3y 

„ 3Zx  3Zy 

z 3x  3y 


•1  2 ‘ 


8Z  3Z 

( — - + — 2-) 
v 3y  3x  * 


This  completes  the  basic  formulation  of  the  laminate  theory. 
As  it  can  be  seen  from  the  equations  Just  developed,  the  re 
suits  involve  numerous  geometric  and  material  parameters. 


I 


i 


A FIVE-LAYER  LAMINATE  WITH  A CRACK  IN  THE  CENTER  LAYER 


Consider  a symmetrically  layered  composite  plate  which 
contains  a crack  of  length  2a  through  the  middle  layer.  The 
five-layer  plate  Is  assumed  to  be  made  of  three  different 
elastic  materials  with  the  stacking  sequence  as  that  shown  In 
Figure  1.  The  subscripts  1,  2 and  3 attached  to  the  Young's 
modulus  E,  Poisson's  ratio  v and  layer  thickness  h refer  to 
those  properties  belonging  to  the  Individual  layer.  In 
Figure  1,  the  crack  length  2a  Is  taken  to  be  small  in  compar- 
ison with  the  in-plane  dimensions  of  the  plate  such  that  the 
boundary  effect  on  the  stress  field  can  be  neglected. 


A remote  tensile  stress  of  magnitude  <jq  applied  to  the 
laminate.  By  the  application  of  the  principle  of  superposi- 
tion, the  equivalent  problem  to  be  solved  here  is  that  of 
normal  traction  applied  to  the  surface  of  the  crack.  This 
traction  can  be  estimated  from  the  same  laminate  without  the 
crack  as  [3]: 


1-V jVm 


c l <Vhl)] 

n-1  n x m-1  al  1 vm 
3 E 2 h 2 l-v* 

< I [(*“■)  (c*)  (rrrrm  + d 

n-1  E1  nl  1 vn 


Ka. 


(13) 


in  which 


E-  h-  1-v, v-  E-  h-  l-v,v- 

D - 2<S7><h7><T^>  * + 


+ 


(l-v12 ) (l-v2v-) 


(14) 


* 


i 


j 


12 


u 


Due  to  symmetry,  the  differential  equations  (9)  are  to  be 
solved  subjected  to  the  following  conditions: 

txyu,°,s>  • TyzU,0,s)  - 0 

o„(x,o,z)  ■ - Ko^fL'(z),  0<x<a 
y * oi  ” 

if  (x,o,z)  » 0,  x>a 

Txy(°*y"l)  ‘ txz(o’y'x>  " 0 * <V°>y'z)  ' 0 
,xy(x,y,o)  ■ Tyz<x'y'o)  ’ °!  Vx>y’o)  • 0 

Eliminating  the  z-dependence  in  equations  (15)  renders 
sxy(x>0)  ’ Vx<°>  * 0 

Sy ( x , o ) ■ - Kp0,  0<x<a  (16) 

u (x,o)  ■ 0,  x>a 
* 

SjjyCo.y)  » Zx(o,y)  - ux(o,y)  - 0 

In  the  absence  of  out-of-plane  bending,  the  last  of  equations 
(15)  is  satisfied  when  the  functions  fp(z)  are  required  to  be 
symmetric  about  the  mid-plane  z"0.  In  equations  (16),  the 
equivalent  stress  pQ  is  introduced  such  that  pQ  ■ Cgh.^2. 


(15) 

I 
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Solution  P/tocedu4.e.  Fourier  transform  techniques  will  be 


used  to  solve  the  system  of  differential  equations  (9). 

Depending  on  the  evenness  and  oddness  of  the  function  with 
respect  to  its  argument,  either  the  Fourier  cosine  or  sine 
transform  may  be  employed.  For  example,  the  function  Sx°(x,y), 

S*  (x,y)  will  denote,  respectively,  the  Fourier  cosine  and 
sine  transform  of  the  functions  S (x,y)  and  S (x,y).  The  same 

a Ay 

notation  applied  to  the  other  functions  S„c(x,y),  Z_S/_  » ...  , 

y A\A^y^j©uc« 

i.e. , 


Sx(x,y)  ■ j / S c(s,y)cos(sx)ds 


S (x,y)  « 7 / S (s,y)cos(xs)ds 
j T o y 


S (x,y)  - 7 / S s(s,y)sin(sx)ds 


xy 


* o xy 


o 00 

ZT(x,y)  * 7 / Zx8(s,y)sin(sx)ds 

o 


Z (x,y ) • 7 / Z c(s,y)cos(sx)ds 


* o y 


(17) 


Z*(x»y)  ■ 7 / Z °(x,y)cos(sx)ds 
o 


and  the  same  applied  to  the  displacements 


ux(x,y)  » 7 / ux8(s,y)sin(sx)ds 
o 

uy(x,y)  « 7 / uy  (s,y)cos(sx)ds 


u,(x,y)  ■ 7 / u c(s,y)cos(sx)da 
o 
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(18) 


i 


■ 


Making  use  of  equations  (17)  and  (18)  and  without  going  into 
details,  a solution  may  be  written  as 

SxC(a»y)  • qB(s)e"qy  ♦ Ijilj  ReCCljtsI^ 

+ pI2P2)  - I4(I1-I2)(sP1-pP2)]e“py} 

“ rpij  + I2p2^Aie”py^ 

SyC(8»y)  " - ^(®)®-qy  - fpij  Re{[I3(sI2P1 

♦Pl^)  + WV (sPx  - pP2)]e’py) 

+ ^zijz  Re[(Ijp2  + I2s2 )A1e“py] 

S%yS(l’y>  - T^rq  XetspAje-^]  - 

+ ReC(pPx  + sP2)e“py]}  (19) 

Zxs(s,y)  ■ I*  B(s)e“qy  + Re[P1e“py] 

Zy8(s,y)  - B(s)e“qy  + Re[P2e"py] 

ZzC(s»y)  ■ ReC(sP1  - pP2)e_py] 

ux8(sfy)  - Re[(2sA1  - I^)*"07]  - I3  f B(s)e“qy 
uyc(a,y)  ■ *91(20*!  - I3P2)e"py]  - I3  B(s)e"qy 
wzc(s,y)  ■ 2 Re[A1e~py] 
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in  which  Re  denotes  the  real  part  of  a complex  function  and 
A^(s)  and  B(s)  are  the  unknowns  to  be  determined  from  the 
boundary  conditions.  In  equations  (19),  the  quantities  p,  q, 
and  P2  represent 


a, 


Cs2  - j--  + i/ 

2a3  a3 


aU 


1/2 


(s2  - 


Yl  ♦ 1t2) 


1/2 


• oR  + iej 

q • (s2  + |j)1/2  (20) 

a,(s2-p2)-a2 

P1  * 28  ag(a'2'-"p2 J+i  Ai 

a1(s2-p2)-a2 

p2  " 2p  i^i^^T+l  " A1 

Enforcing  the  boundary  conditions  (16)  yields 

B(,)  * - ReC2o  *i(»n 

(21) 

Im[Ai(s)]  - 0(8)  ReCA1(s)] 


as  well  as  a set  of  dual  integral  equations 

2*i 

r It  A(s)cos(sx)ds  * 0,  x>a 
* o 8 

(22) 

2 * 

r / g(s)  A(s)cos(sx)ds  * - Kp  0<x<a 
o o - 

governing  the  unknown  A(s)  -where 

A(s)  • sg-^s)  Re[A1(s)]  (23) 
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I 


1 


In  equations  (21),  (22)  and  (23)  the  following  contractions 
have  been  made: 

0(s)  - - [pR(s2  + ^ «i)  + Pj  ^ 623/CpR  2^  ^ 


' pI(s2  + 2a£  61)] 

^alYl“a2^a6Yl+l3  + ala6Y22 
(a6Yl+^2  + a**Y92 


6 y2 


. T2(yaaV 

(a6Yi+1J*  + 


gx(s) 


g2<3) 


2 Re[p (l+i0 ) ] 

ipij  R.((l*i»>CCI1P*«2.»>  * 


(24) 


__2 

+1 


tl3(I1+I2)qp  - 


I3(I2s2+  Ijp2) 


- i4(i1-i2)(s2-p2)]1]} 


g2(s) 

s(s)  ■ SS^T 

Following  [1],  the  solution  to  equations  (22)  may  be  expressed 
in  terms  of  a new  function  «(e)  as 
KP«a  i 


A(s)  - — £ — s / 4%  *(C)  J0(sac)dc 


(25) 


where 


(1---j0  3)(I1“I2) 


{(Ij+  I2)  + 63  C-1^2  + <Ii“I2)^  " I4) 


(▼l«2+Y2«i)  “ «2  Sa£  (I1+I2)]} 
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(26) 


sxy(*,y> 


2XP''a  "t^To  Reli^§7  C(iti8)<i-(s1-i«2)i3»e 


/ 


-py 


I,(s2+q2) 

+ ( 6 1— ± <5 2 ) ( H-i 8 ) g3z’ — e“qy ] }sin(sx)ds 


(30) 


Zx(x,y) 


2Kp  a2  - 


/ Re{(61-i52)(l+i6)(se"py 


- f^e-^)  }sin(sx)ds 

9 

2Kp  a2  - . 

Zy(x,y)  - -K°-  ■ f Re{(«1-I«2)(l+ie)p(e"py-e“qy)} 


cos (sx)ds 

2Kp  a2  - ... 

Zz(x»y)  - — K9-  ! Re{(«1-H2)(l+iB)(r1-iY2> 


e”py } 


cos (sx)ds 


In  which 

I 

*(»)  - / /€  *(€)J0(sa5)d^  (31) 

Equations  (30)  in  conjunction  with  the  thickness  functions 
fp(z)  determine  the  stress  distribution  in  the  laminated  plate. 

Thin- tUme.ru iona.1  Cxack  F nont  S& te**  field.  The  singular 

behavior  of  the  quantities  S , S , etc.,  can  be  obtained  by 

* y 

application  of  equations  (30)  in  which  only  the  terms  that 
contribute  to  the  divergence  of  the  improper  integrals  at 
x-±a  and  y*0  are  retained.  Expanding  the  integrals  in 
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equations  (30)  for  large  s and  multiplying  the  stresses  by  the 
functions  f-^z),  f[(z)  and  f![(z)  according  to  equations  (1), 
the  following  crack  front  stress  distribution  is  obtained: 

o ( z)a1/2  , 

‘ K,(1)  ? ‘ 1/2  [oos  2 - 3 f-3  + 0(r°) 


a (z)a1/2  - I-,- 

°y  ' K‘(1>  ' 1/2  [co°  | ft]  ♦ 0(r«) 


t * K*(l)  — -■  i ■ i sine  cos  4—  + o(r°) 

xy  (2r)1/2  2 2 


- K«(l) 


aT(z)a1/2 

' 77,.  a/2  cos  2 + 0(r  > 


t * t * o(r  ) 
xz  yz  v ' 


In  equations  (32),  r and  e are  the  local  polar  coordinates  with 
origin  at  the  crack  point  x-a  and  y«o  in  the  xy -plane, 
symbol  *(1)  is  the  value  of  the  function  ♦(?)  as  *+l,  the 
location  of  the  normalized  crack  tip.  The  functions  <jp(z)  and 


aT(z)  are  defined  by  the  expressions 

c w - jf  q (z.)p0 


,(z)  f ^( z )p( 


in  which 
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M 


( 1“2^6 3)  1 [(I1+I2}  2^  (5562+665l}  + T^l  ~T&£ 


N - (l-^^r1  [53(6251+6162)] 


(34) 


The  customary  inverse  square  root  singularity  at  the  crack 
front  and  the  singular  variations  of  the  stresses  are  observed 
in  equations  (32).  This  is  in  agreement  with  the  exact  three- 
dimensional  crack  edge  solution  [4].  The  approximation  of 
this  theory  enters  only  through  the  intensity  of  the  local  stress 
field.  In  this  connection,  the  stress  intensity  factor  may 
be  defined  by 


k. (z)  ■ K«(l)  a (z)a 
1 P 


1/2 


(35) 


* 


for  the  in-plane  stresses  and 

kx(z)  - K*(l)  oT(z)a1/2  (36) 

for  the  transverse  stress  az.  Up  to  this  point,  the  in-plane 
behaviors  of  the  stresses  are  completely  determined.  Additional 
conditions  based  on  the  characteristics  of  the  stresses  in  the 
interior  region  of  the  layer  and  near  the  free  surface  on 
interface  are  required  for  the  determination  of  the  functions 
f^(z).  Photoelastic  Analysis  [4]  has  shown  that  the  stresses 
in  the  interior  of  plate  are  nearly  in  a state  of  plane  strain 
while  they  undergo  large  variations  in  a layer  near  the  plate 
surface.  Substituting  equations  (32)  and  (33)  into  the  plane 
strain  condition  yields 

f1(z)  - A cos  (p2z)  _2i-  (37) 


The  parameter  controls  the  qualitative  feature  of  the  stress 
distribution  across  the  layer  thickness  and  has  the  form 

N 1/2 

Pl  - <2^m>  (38) 

The  constant  A in  equation  (37)  serves  as  a quantitative 
measure  of  the  stresses. 


Mo  ruing  ulLolk  S£> tea*  Vl4£/Ubu£lon  A c.ko&/>  The.  Laminate..  It  is 

observed  from  equations  (5)  that  the  stresses  in  the  laminate 

can  be  found  by  multiplying  the  appropriate  functions  fp(z)  or 

its  derivatives  with  the  in-plane  quantities  S , S etc., 

X y, 

defined  in  equations  (30).  Since  S , S etc.  may  be  evaluated 

x y j 

by  numerical  integration  and  the  numerical  solution  of  the 
Fredholm  integral  equation  (38),  the  task  here  is  to  construct 
the  functions  fg(z)  and  f^(z)  that  are  consistent  with 
equations  (3)  and  (37).  To  this  end,  the  following  expressions 
for  fp(z)  may  be  written: 

f1(z)  ■ A cos  (p1z),  O^z^i^ 


pi 

f2(z)  - A{cos  (p-jh^  cos  [p2(z-h1)]  - — sin  (p1h1) 


p2 


f3(z) 


sin  [p2(z  -h1)]},  h^z^h^hg 
A{ [cos 

(P2h2)] 


(39) 


pl 

:os  (p1h1)  cos  (p2h2)  - ~ sin  (p^h^  sin 

i2)]  cosCp^z-t^-hg)  ] - ~ [p2  cos  (Pjh^  sin 

- '•  )•  **  : ■ 

u)  + P2  sin  (p^)  cos  (p2h2)]  sin  [p^ 


-Vi 


V ^ . 


u 


( 

! 

I 


The  parameters  Pj » J - 1,2,3,  are  defined  as 


T 


(40) 


With  these  definitions  for  fp(z),  the  stress  field  in  an 
arbitrary  point  in  the  laminate  may  be  obtained  from  the 
following  expressions : 

«x<X>y,a)  . f"  (Z)  / lifly  R.(Il82+I2p2) 

- l]e“py  - (61-i62)(l+i6)(Y1-iY2) 

1 +X  e °y  “ ^ 5, -i62 ) ( 1+ig ) j — 3— e”qy)cos  (sx)ds 
12  l*x2 

2 Kp  fit  ^ , 

oy(x,y,z)  - -f-  r;<z)  / iifiy  Re{ (I2s2+I^p  2 ) ^f|I 

[l-(6^-i62 Jl^le  py  - (6^-152 ) (1+iB ) ( Yx“lY2 ) 

1 K __  1 <3 

T +t  e P + (ii-lfip) (1+iB)  y — f-e_qy}cos  (sx)ds 
12  c j‘1~‘l2 

2Kp  a2  - , . 

Txy(x,y,z)  - — ^ — f"(z)  / Re<T^§lj  Cd+iB) 

rtU  I,(s2+q2) 

<l-(61-U2)I3>.e  py  +( 61-i«2 ) (1+iB)  -3  ^s2 

e^^jjsin  (sx)ds 


-23- 


Because  of  assumptions  (2,),  (3)  and  (4)  the  stress 
expressions  (41)  are  not  expected  to  give  satisfactory  results 
In  the  regions  bounded  by  x»a,  y»o  throughout  the  laminate 
except  on  the  middle  layer  that  contains  the  crack.  Moreover, 
as  mentioned  In  [1]  and  [5]  the  plane  strain  condition,  which 
leads  to  the  construction  of  the  functions  fp(zj,  are  not 
expected  to  hold  near  a free  surface  or  Interface  because  of 
the  lack  of  mechanical  constraint.  Hence,  fp(z)  must  be 
constructed  separately  whenever  the  condition  «jz»  v(ax+oy)  is 
not  satisfied.  The  concept  of  a boundary  layer  can  then  be 
Introduced.  Such  a procedure  results  In  additional  discussions 
and  Is  beyond  the  scope  of  the  present  Investigation. 
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NUMERICAL  RESULTS  AND  DISCUSSIONS 


The  numerical  results  in  this  report  are  divided  into  two 
parts.  The  first  pertains  to  the  five-layer  laminated  compo- 
site where  attention  is  focused  on  the  variations  of  the  stress 
intensity  factor  with  the  material  geometrical  parameters.  The 
second  considers  the  stress  distribution  in  a laminated  plate 
due  to  the  influence  of  adhesive  layers.  This  is  done  for  a 
specific  example  of  aluminum-epoxy-steel  system.  The  stress 
intensity  factor  as  well  as  stress  field  away  from  the  crack 
are  determined  for  this  example.  For  convenience,  the  Poisson's 
ratios  for  all  the  numerical  computations  are  taken  to  be 
0.3. 


The.  F-tve-  Lnyex  PxobZem.  The  stacking  sequence  of  a laminate 
is  determined  by  assigning  specific  values  of  the  material  and 
geometrical  parameters  and  will  affect  the  intensity  of  the 
stresses  transmitted  to  the  crack  front  region.  The  stress 
intensity  factor  gave  a measure  of  this  influence  and  in  its 
normalized  form  is  given  by 


kx(z) 


17? 


op(z)a 


K#(l) 


(42) 


in  which  K is  given  in  equation  (13)  and  «(1)  may  be  calculated 
from  the  Fredholm  integral  equation  (28).  The  variation  of  1E1  as 
a function  of  the  center  layer  thickness  to  crack  length  ratio 
hj/a,  is  plotted  in  Fiqure  2 for  various  values  of  P ■ p^h^. 

The  Young's  modulus  and  thickness  ratio  are  taken  to  be 


* 
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E2/Ex  ■ 0.1,  E^/E^  ■ 1.0,  h2/h1  » 0.1  and  h^/i^  ■ 1.0.  The 
choice  of  P depends  on  the  transverse  stress  distribution  in  a 
given  problem.  Numerically,  it  is  confined  to  0<P<|-  because 
of  the  cos  (P)  dependence  of  the  stress  field  which  changes 
sign  when  P is  outside  of  the  above  range.  Note  that  as  P is 
decreased,  k^  rises  more  sharply  for  small  values  h^/a.  This 
type  of  behavior  was  also  observed  for  the  three-layer  problem 
Cl]. 


I 

r , 

r- 


Pigure  3 shows  a plot  of  versus  h^/a  for  “ 0,1>  1,0 

and  10.0.  The  parameter  P is  chosen  to  be  0.7  while  ^2^1  “ 0,1» 
h2/hi  ■ 0.1  and  h^/l^  » 5.0.  The  stress  intensity  factor 
increases  with  decreasing  stiffness  ratio  E^/E^.  The  same 
trend  is  observed  in  Figure  4 and  5 where  the  thickness  ratio 
h^/h^  is  taken  to  be  1.0  and  0.2,  respectively.  By  comparing 
the  above  three  figures,  it  is  evident  that  the  stress  intensity 
factor  increases  with  increasing  h^/h1  for  E^/E1<1.0.  On  the 
other  hand,  the  opposite  trend  is  observed  for  E^/E^l.O  This 
behavior  is  more  clearly  illustrated  in  Figure  6 where  the 
variations  of  k^  versus  hg/i^  for  three  E3/E 1 ratio  are  plotted. 
The  values  of  the  other  parameters  are  fixed  at  P ■ 0.7, 

Eg/^i  ■ 0.1,  h2/h^  ■ 0.1  and  h^/a  * 1.0.  Thus  the  crack  tip 
stress  can  be  lowered  by  (a)  raising  the  stiffness  ratio  £3/2!* 
or  (b)  decreasing  the  thickness  ratio  h^/l^  for  E^/E^l.O,  or 
(c)  increasing  the  h^/h^  ratio  for  E^/E^l.O.  These  conclusions 
are  valid  as  long  as  the  parameters  in  layer  number  two  are 
kept  constant. 
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The  variations  of  kx  with  h1/a  are  displayed  In  Figure  7 
for  fixed  values  P ■ 0.7,  " 1,0 > h3/hi  “ 1*0  and 

hj/hi  ■ 5.0.  The  three  curves  represent  stiffness  ratios  of 
E2/E1  ■ 0.1,  1.0  and  10.0.  Hence,  ic^  decreases  with  Increasing 
E2/Ei  ratio.  Similar  trends  are  illustrated  in  Figure  8 
for  h2/h1  ■ 1.0  and  in  Figure  9 for  h2/h1  » 0.1.  The  effect 

of  the  thickness  ratio  hj/h^  on  the  stress  intensity  factor 
is  demonstrated  in  Figure  10  for  P • 0.7,  E^/E 1 ■ 1.0,  h^/l^  » 
1.0,  h1/a  ■ 1.0.  Therefore,  the  same  conclusion  as  that  for 
the  outside  layer  can  be  drawn  here.  The  stress  intensity 
factor  decreases  with  (a)  increasing  E2/E1  ratio,  or  (b) 
decreasing  hg/l^  ratio  if  E2/E1<1.0,  or  (c)  increasing  h2/h1 
ratio  if  E2/E1>1.0. 

Cxacked  Laminate.  With  Adhesive  Layext . As  mentioned 
before,  the  five-layer  plate  problem  may  be  used  to  model  the 
presence  of  adhesive  layers  in  laminated  structures . Let  mate- 
rial 2 in  Figure  1 represent  the  adhesive  layers.  As  a specific 
example,  numerical  computations  were  performed  on  an  alumlnum- 
epoxy-steel  system.  Layer  one  is  assumed  to  be  made  of 
aluminum,  layer  two  of  epoxy  and  layer  3 of  steel.  Therefore, 
the  stiffness  ratios  are  E^/E^B.O  and  Eg/E^^  • 0.05.  The 
thickness  ratio  h^/h^  is  taken  to  be  one.  Since  the  adhesive 
layer  is  about  one  order  of  magnitude  smaller  than  the  other 
layers,  the  thickness  ratio  h2/h1  should  be  such  that 
0<h2/h1<>0. 1.  Within  this  range,  the  scales  of  Pigure  10  show 


that  h2/h1  has  little  effect  on  the  stress  intensity  factor. 
This  is  shown  by  the  ^ versus  h2/h ^ in  Figure  11  for  the 
aluminum-epoxy-steel  system.  The  crack  is  in  the  aluminum 
layer.  The  parameters  are  fixed  at  P ■ 0.7,  h^a  » 1.0  and 
v1*v2  ■ v ^ * 0.3* 

It  should  be  said  that  although  the  adhesive  layers  have 
negligible  influence  in  the  stress  intensity  factor,  their 
presence  may  affect  inter  layer  separation  or  delmaination. 

For  this  reason,  it  is  essential  to  turn  to  the  calculation  of 
stress  distribution  in  the  adhesive  layers. 

Define  the  normalized  stress  7 ■ a/o„  where  o is  the 

o o 

applied  stress.  Figure  12  depicts  the  variations  of  the 

stresses  as  with  the  normalized  thickness  parameter  z/h ^ for 

h2/hi  ■ 0.0  or  for  a three-layer  laminate.  The  stresses  are 

computed  at  the  point  x/a  - 1.0  and  y/a  -3.0.  It  is  seen 

that  7X,  Oy,  az  and  7xy  decrease  in  magnitude  as  z/h ^ is 

increased.  The  maximum  value  of  these  stresses  occur  at  the 

center  of  the  plate,  l.e.,  z ■ 0.  For  the  transverse  shear 

stresses  7„„  and  7.  . the  trend  is  reversed.  A similar  plot  is 
xz  yz 

Figure  13  for  h2/h1  » 0.1  exhibits  the  same  behavior  except  the 
magnitudes  are  slightly  smaller  than  their  counterpart  in 
Figure  12.  In  Figure  14,  variations  of  the  stresses  for  h2/h1 

■ 1.0  are  displayed.  Unlike  the  previous  two  cases  for  h2/h1  ■ 
0.0  and  0.1,  cv,  o„,  7_„  and  a changes  sign  when  the  free 
surface  is  approached.  Figure  15  depicts  the  stress  variation 
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across  the  thickness  at  a different  point  x/a  ■ 3.0  and  y/a  ■ 

3.0.  This  point  is  further  away  from  the  crack  as  compared 

with  the  point  in  the  previous  figures.  For  h2/h1  ■ 0.0,  the 

stress  components  0 a , 7 and  7 decrease  in  magnitude 

x y xy  z 

as  z/h,  increases.  However,  o becomes  compressive.  7__  and 
1 y xz 

7yz  again  rises  with  increasing  z/h1  except  has  negative 
sense  while  7yz  is  positive.  Similar  plot  for  h2/h1  - 0.1  is 
shown  in  Figure  16  for  comparison.  Figure  17  displays  the 
variation  of  7j  versus  z/h1  for  hj/h^  ■ 1.0.  It  is  interesting 
to  note  that  the  transverse  shear  stresses  reach  a peak  and 
then  drop  off.  Figure  12-17  give  some  insight  into  the  stress 
distribution  across  the  thickness  in  the  laminate. 

To  illustrate  the  stress  variation  in  the  xy-plane  at  a 

given  z elevation,  three-dimensional  plots  of  the  normal  stresses 

7_,  a„  and  a are  depicted,  respectively,  in  Figures  18,  19  and 
x y z 

20  as  functions  of  the  position  parameters  x/a  and  y/a  for 
z/h^  ■ 1.  These  stresses  are  thus  corresponding  to  the  inter- 
facial stresses  between  layer  one  and  three  since  b2/h1  ■ 0 in 
these  figures.  It  is  observed  that  all  the  stresses  die  out 
quickly  with  increasing  x/a  and  y/a  values.  All  three  stress 
components  changes  sign  from  positive  to  negative  before 
reaching  zero  asymptotically. 

More  detailed  stress  variations  as  a function  of  x/a,  y/a, 
z/h,  and  h2/h^  are  tabulated  in  Tables  1-9 • Table  1 shows  the 
stress  field  for  a three-layer  laminate  at  the  middle  of  the 


center  layer.  The  parameters  x/a  and  y/a  are  extended  from 
1. 0-5.0  and  1.0-10.0,  respectively.  No  need  to  expand 

these  ranges  since  the  stresses  have  already  reached  their  , 

asymptotic  values.  Table  2 exhibits  the  stresses  at  i 

the  Interface.  The  stress  field  In  the  middle  of  the  outer 
layer  Is  listed  In  Table  3.  By  taking  h2/h1  ■ 0.02,  Table  4-6 
shows  the  stress  distribution  as  a function  of  x/a  and  y/a  at  the 

i 

middle  of  the  center  layer,  of  the  adhesive  layer  and  of  the 
outside  layer.  Finally,  Table  7-9  lists  the  variations  of  the 

i 

stresses  for  h^/h^  ■ 0.1  at  the  middle  of  each  layer. 

To  recapitulate,  the  stress  calculations  present  here  are 
essential  to  the  prediction  of  failures  by  the  delamination 
mechanism.  A failure  criterion  must  be  used  in  conjunction 
with  these  stress  results  to  forecast  prediction.  This  phase 
of  the  work  will  be  left  for  future  research. 

! 

] 


* 
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CONCLUSIONS 

A method  for  computing  the  stress  distributions  in  a five- 
layered composite  plate  was  developed.  By  knowing  the  fracture 
toughness,  k^Q,  of  the  middle  layer  which  contains  the  crack  the 
load  at  instability  may  be  estimated.  To  check  for  delamination, 
the  stresses  obtained  in  the  adhesive  layers  may  be  used  to  compute 
for  the  strain  energy  density  factor,  S [6].  The  location  of  mini- 
mum S or  Smin  predicts  where  failure  may  occur  while  the  critical 
value  of  Scr,  indicates  when  failure  may  occur.  The  critical 

load  obtained  in  this  way  may  then  be  compared  with  that  found 
from  the  critical  stress  intensity  factor  calculation  to  determine 
whether  or  not  de lamination  has  occurred  before  crack  propagation. 

On  basis  of  the  aluminum-epoxy-steel  laminate  example, 
several  points  of  interest  may  be  concluded: 

(1)  In  the  case  of  a five-layered  laminate,  the  stiffness 
ratios  E2/E1  and  ® * and/°r  the  thickness  ratios  h2/h^  and 
h^/l^  may  be  combined  in  such  a way  to  lower  the  crack  front 
stress  intensity.  This  can  elevate  the  allowable  design  stress 
and/or  prolong  the  life  expectancy  of  the  composite  structure. 

(2)  The  adhesive  layers  in  the  laminate  do  not  appear  to 
have  a significant  influence  on  the  stress  intensity  factor  for 
adhesives  in  the  practical  thickness  ratio  range  of  0<h2/h^<0.1. 
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(3)  The  present  scheme  of  stress  analysis  coupled  with 
a realistic  fracture  criterion  may  be  stored  In  the  computer 
for  otlmizing  the  design  of  laminates  against  brittle  fracture. 
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Figure  3.  - k,  vs.  h,/a  for  various  ratios 
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Figure  9.  - k,  vs.  h,/a  for  various  ratios  of  Ep/E,  (hp/h 
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Figure  17.  - Variation  of  stresses  across  th  thickness 
for  h?/h.  - 1.0  (x/a  - 3.0,  y/a  3.0) 
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Figure  20.  - Variation  of  a as  a function  of  x/a  and  y/a 


TABLE  1 



p = 

. 700  NU 

= .3  CO  Hl/A 

= .500 

E3/E 

Is  3.000 

H3/H1=  1.00 

E2/E1=  . 

050  H2/H1= 

0.00 

Z/Mi 

= .50  _ 

X/A  = 1.00 

- - 

— 

- - — — 

--  — 

y/a 

sx 

SY 

TXY.  _ 

• 

. sz  . 

TXZ 

TYZ  ...  ... 

1.0 

6.77399 

9.05259 

9.93543 

-9.91788 

-3.37776 

•21446 

2.0 

2.  16804 

4.41560 

2.36339 

-4.02311 

-.91718 

-.12796 

3.0 

.79  579 
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.62581 
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-.24877. 
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4.0 

. 30502 

.66399 

.18164 
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-.07292 

- .02395 

5.3 

.1185  8 
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.05621 

-.228  01 

-.02252 
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6.3 
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.01782 
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7.0 
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• .03770 

.00617 
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-.00250 

- .00107 
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.01439 

.00206 

-.01399 

-.00371 

- .00036 

9.0 

.00274 

.00548 

.00069 

-.0  05  60 

-.00036 

0 .00000 

10.0 

.03069 

.00236 

0.00000 

-.0  0140 

0.00000 

0 .00003 

Z/Hl 

= .50 

X/A  = 2.00 

y/a 

sx 

SY 

TXY 

SZ 

TXZ 

T >Z 

1.0 

3.12971 

-1  .93636 

1.50.317 
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_ .1 ..  04371_ 

2.0 
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4.0 
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6.0 
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.00036 
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0 .00000 
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0 .00003 

9.0 
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10.0 
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.00069 
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TABLE  1 (Continued) 
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Z/Hl  * 

.50  X/A 

= 3.03 

Y/A 

SX 

SY 

TXY 

sz 

TXZ 

T >Z 

1.0 

.27623 
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.19444 
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2.0 
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SY 
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SZ 
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2.0 
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3.0 
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P = .700  NU  = .3  CO  ~H1/A  = .500 


E3/E 

1=  3.000  H3/H 1*  1.00 

E2/E1*  . 

050  H2/H1= 

0.00 

7/HI 

= 1.00  X/A 

= 1.00 

• 

- - . 

- 



• 

y/A 

sx 
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sz 

TXZ_ 

TYZ  . 

1.0 
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SX 

SY 
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sz 
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p = . 

700  MJ  ■ 

* .3  00  Hl/A 

* • 500  ” 

E3/E1 

* 3.000 

H3/H1*  1.0  0 

E2/E1*  • 

050  HZ/ HI* 

0.00 

I/Hi 
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* ■ - — 

— — 

Y/A 

SX 

SY 

TXY  _ . 

sz  . ... 
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_T_>Z  

1.0 
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2.0 
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1.25185 
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-2.32017 
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3.0 
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-.820.97  . 
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_ - .16095  _ 

4.9 
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• .06058 

5.0 
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.13506 

.02977 
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-.05696 

- .02170 

6.0 

. 02469 

.05192 

.00944  . 

*.  -.04668  . 

. -.01899. 

-.00723  _ _ 

7.0 

.00983 

.01997 

.00327 

-.01778 

-.00633 

- .30271 

8.9 

. 0 0 363 

.00762 

.09109 

-.007  41 

-.00131 

- .30090 

9.0 

. 0 0 145 

.00290 

.90036 

-.00296 

-.00090  . 

0 .00009 

10.3 

. 00036 

.00109 

9.00000 

-.00074 
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0 .00000 

Z/Hl 
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X/A  = 2.30 

V/A 

SX 

SY 

TXY 

SZ 

TXZ 

T >Z 

1.0 

1.65776 

.-1 .02566 

..  .79620. 

-.67352 .. . 

_-JL.  83914. 

_ 2 .,640 26  

2.0 
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.65  824 

-.65871 

-1.32917 

.74054 

3.0 
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.25269 

.26940 

-.35714 

-.53376 

.16095 

4.0 

.11618 

.14341 

.09766 

-.16449 

-.19169.  . 

.03255 

5.0 

. 04611 

.06717 

.03449 

-.07113 

-.06781 

.00633 

6.0 

.01852 

.02941 

.01234 

-.02519 

-.02351 

.00090 

7.0 

. 03726 

.01234 

. • 00436„ 

. -.01260 

-.00814  . 

. 0.00000 

S.O 

.00293 

.00436 

• 00 145 

-.00519 
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0.00000 
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. 00109 
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•00073 
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-.00090 

0.00030 
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.00036  . 
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TABLE  3 (Continued) 
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I 

- , 

Z/Hl 

= 1.50  X/A 

= 3.00 

.1 

• 
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sx 

SY 

TXY 

sz 

TXZ 

~ Tyz  j 
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3.0 
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4.0 
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5.0 
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6.0 

. 00  472 

•00109 

.00363 
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.00452 

7.0 

. 03218 

.00109 

•00182 
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8.0 
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.00090 
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. 00036 
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-.00074 

, -.00390 

, J).. 00000.  - - . 

10.0 

.00036 

.00036 

0.00000 

0.00000 

0.00000 

0 .00000 

Z/Hl 

= 1.50  X/A 

= 5.00 

Y/fl 

sx 

SY 

TXY 

SZ 

TXZ 

T >Z 

1.0 

-.  02578 

.00508 

-.01307 

-.00741 

,.05968 

-02170  

2.0 

-.01525 

-.01271 

-.01271 

.00889 

.03707 

- .00271 

3.3 

-.  00363 

-.01380 

-.00508 

.00667 

.01356 

.00814 

4.0 

.00109 

-.30799 

-.00036, 

•00222 

.00181. 

. .30814 

5.3 

. 0 0 18  2 

-.00327 

.00109 

0.0  00  00 

-.00181 

.00452 

6.0 

.03109 

-.00109 

.00073 

-.00074 

-.00181 

.00271 

7.0 

. 00073 

-.00036 

.00073, 

-.00074 

,-.00090 

.30090 

8.0 

.03036 

0 .90000 

.00036 

0.00000 

-.00090 

0 .00000 

9.3 

.00036 

0 .00300 

0.00000 

0.000  00 

0.00000 

0 .00000 

10.0 

0.03003 

0 .00000 

0.00000 

_ 0.00000 

0.00030 

.0  .00000 

-59- 


T 


TABLE  4 


P s .700  NU  * .3  CO  Hi/ A * .500 


E3/E1 

= 3.000 

H3/H 1*  1.00 

E2/E1*  . 

050  H2/H1* 

.02 

Z/Ml 

= .50 

X/A  *..1.00 

• 

_ _ . 

■ - . — _ _ 

Y/A 

sx 

SY 

_ TXY  . 

sz_._ 

_TXZ  ...  . _ 

1.3 

6.84676 

9.14311 

10.03479 

-10.01706 

-3.41154 

2.0 

2.18972 

4.45975 

2.38703 

-4.06334 

-.92635- 

3.3 

.80375 

1.76742 

.63206 

..-1*  5 65.43 

-.25126 

4.0 

.30807 

.67568 

.16346 

-.599  05 

-.07365 

.5.3 

.11977 

.25753 

.05677 

-.23029 

-.02274 

6.3 

. 0 4 70S 

.09900 

.01800 

-.08901 

-.00758 

7.3 

. 01869 

.33608 

.00623 

-.03391 

-.00253 

8.3 

.00  692 

.01454 

• CO  208 

-.01413 

-.00072 

9.3 

.00277 

.00554 

..  • 0C069 

-.00565 

-.00036 

10.3 

i 

.00069 

.00208 

0.03000 

-.00141 

0.00000 

• •- 

Z/Hl 

* .50 

X/A  * 2.00 

Y/A 

SX 

SY 

TXY 

SZ 

TXZ 

1.3  . 

3.16101 

-1  .95573 

_ 1.51820 

-.1.28427  _ 

_-.73429 

2.0 

1.46420 

.33507 

1.25513 

-1.25602 

-.53066 

3.3 

.4  9 291 

.48164 

.51366 

-.68099 

-.21191 

4.0 

.22153 

.27346 

.18623 

-.3.1365 

-.0765.3 

5.0 

. 08792 

.12807 

.06577 

-.13563 

-.02708 

6.0 

. 03531 

.05608 

.02354 

-.046  04 

-.00939 

7.3 

. 01385 

.02354 

.00831 

-•024C2 

-.00325 

8.3 

. 00554 

.00831 

.00277 

-.00989 

-.00108 

9.3 

.00208 

•0C415 

.00138 

-.0  0424 

-.00036 

10.0 

. 00  069 

.00138 

.00069 

-.00141 

0.00000 

I 

II 


_ T YZ 

.21661 

- .12924 

- .06426 

- .0  2419 
• .00666 

.-..00289 

- .00108 
- .00036 
0 .00000 
0 .00000 


T >Z 

1.05415  ..  ... 
.29567 
.06426 
_ .0130  0 
.00253 
.00036 

0 .00000  . . . 
0 .00000 
0 .00000 
0 .00000 
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TABLE  4 (Continued) 


Z/Hl  a .50  X/A  a 3.00 


Y/A 

sx 

SY 

TXY 

sz 

TXZ 

T >Z 

1.0 

.27899 

-.81621 

-.02492 

.19639 

-.00650 

.20325 

2.0 

.39115 

-.35930 

•21461_ 

-.10596 

-.10325  _ 

.15199 

3.0 

.22776 

-.04500 

.17515 

-.148  35 

-.07581 

.0  6245 

4.0 

.10  eoo 

.03392 

.09069 

-.10031 

-.03827 

.02130 

5.0 

.04777 

.03254 

•04015 

-.05369 

-.01661 _ 

.00686 

6.0 

.02077 

•01938 

•01662 

-.02543 

-.00686 

.00217 

7.0 

.00831 

.00969 

.00623 

-.01130 

-.00253 

.00072 

8.0 

. 0 0 246 

•00485 

.00277 

-.00565 

-.00108 

..00036 

9.0 

.00138 

.00208 

.00069 

-.00283 

-.00036 

0 .00000 

10.0 

.00069 

.00069 

.00069 

-.00141 

0.GQ0Q0 

0 .00000 

2/HI 

a .50 

X/A  a 4.00 

Y/A 

sx 

SY 

TXY 

SZ 

TXZ_  .... 

_T  >Z 

1.0 

-. 12115 

-.16615 

-.08515 

.09607 

.05704 

00975 

2.0 

.00831 

-.16061 

-.02354 

.05651 

.01264 

.03105 

3.0 

.04  Cl  5 

-.07892 

.01869 

.00565 

-.00758 

.32310 

4.0 

. 33046 

-.02423 

.02146 

-.01272 

- • 0 C9Q  3 

.01155 

5.0 

. 01711 

-.00346 

■ C 1385 

-.01130 

-.00578 

.00505 

6.0 

. 0 0 50  0 

.00208 

.00692 

-.00848 

-.00289 

.00161 

7.0 

.00415 

.00208 

.03346 

-.004  24 

-.00144 

.00072 

8.0 

.00  20  8 

.00138 

.00138 

' -.00283 

-.00072 

.00036 

9.0  . 

. 00069 

„• 00069_ 

•00069 

. -.0  0141 

-.0  0036 

_ 0 .00000 

10.0 

.00069 

.00069 

0.00000 

0.000  00 

0.00000 

0 .00000 

Z/Hl 

a .50  X/A  a 5.0C 

V/A 

. SX 

SY 

TXY 

SZ 

TXZ 

T >Z 

1.0 

-.04915 

.00969 

-.02492 

-.01413 

.02383 

.-  .00866 

?.Q 

-.02908 

-.02423 

-.02423 

.01695 

•01480 

- .00108 

J.O 

-.  00  €92 

-.02631 

-.00969 

.01272 

.00542 

.00325 

t.O 

.00208 

-.01523 

-• 0C069 

.00424 

. .0007.2 

...  .00325 

5.0 

.00  346 

-.00623 

.00208 

0.00000 

-.00072 

.30181 

e.o 

. 00208 

-.00208 

•00138 

-.00141 

-.00072 

.00108 

7.0 

.00  138 

-.00069 

.00138 

-.00141 

->.00036 

...  .33036 

8.0 

.00069 

0 .00030 

•03069 

0.00000 

-.00036 

0 .30000 

9.0 

. 00  069 

0 .00000 

0.00000 

0.00000 

0.00300 

0 .00000 

10.0 

0.00030 

0 .00000 

0.00000 

0.00000 

0.00000. 

0 .00000 
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TABLE  5 


p = . 

700  NU 

= .3  00  Hl/A 

= .500 

• 

E3/E1 

s 3.0  CO 

H3/H1*  1.00 

E2/E1*  . 

050  H2/HI* 

.02 

• 

_Z/H1 

II 

• 

o 

X/A  * 1.00._ 

--  • 

— - - 

— 

- 

Y/A 

sx 

. SY  . 

TXY 

. _sz 

..TXZ 

_T  YZ  ...  . 

1.3 

5.5416ft 

7.40030 

9.12200 

-8.10765 

-6.46253 

.41032 

2.3 

1.77233 

3.60966 

1.93202 

-3.28880 

-1.75490 

- .24482 

3.1 

.65054 

1.43053  . 

.51158. 

-1.26704. 

„ -.47597. 

-.12173 

4.3 

.24935 

.54688 

.14849 

-.48486 

-.13951 

- .04582 

5.1 

. 09  694 

.20844 

.04595 

-.18640 

-.04308 

- .01641 

6.3 

. 03910 

.08013 

.01457. 

..  r. 07204  . 

“_»01436 

.^..00547  . 

7.3 

.01513 

.03082 

.00504 

-.02744 

-.00479 

- .00205 

ft. 3 

.03560 

.01177 

.30168 

-.01144 

-.00137 

- .00068 

9.3 

.03224 

.0044ft 

.30056 

-.00457 

_-.  00068 

0 .00000 

10.3 

. 00  056 

•00168 

0.00000 

-.00114 

0.003-13 

0 .00000 

Z/Hl 

* 1.01 

X/A  = 2.30 

Y/A 

sx 

sv 

TXY 

sz 

TXZ 

T YZ 

1.3 

. 2. 55»47 

.-1.58294 

. 1 . 22881 

-1.03947 

..-1.39098 

1 .99689 

2.0 

1.19510 

.27120 

1.01588 

-1.01660 

-1.00528 

.56009 

3.3 

.39  896 

.38999 

.41577 

-.55118 

-.40143 

.12173 

4.0 

. 17931 

.22133 

.15073 

-.25386 

-.14498 

.02462., 

5.0 

.07  116 

.10366 

.05323 

-.10978 

-.05129 

.00479 

6.0 

. 02958 

•04539 

.01905 

-.03888 

-.01778 

.00068 

7.0 

.01  121 

.01905 

.00672 

-.019  44 

-.00615 

. 0 .00000  _ 

6.0 

.00449 

.00672 

.00224 

-.30800 

-.00205 

0.03000 

9.3 

. 00 168 

.00336 

.00112 

-.00343 

-.00068 

0 .00000 

10.3 

.00056 

.00112 

. .00056 

-.00114 

0.00000 

0.00000 

TABLE  5 (Continued) 


T 


T 


r 


Z/Hl  = 

1.01  X/A 

= 3.00 

Y/A 

sx 

SY 

TXY 

SZ 

TXZ 

T YZ 

1.1 

.22581 

.66063 

-.02017 

.158  95 

-.01231 

.38502 

2.0  . 

.31659 

.29081 

....  17370  _ 

-.08577.. 

-.19559 

28791 

3.0 

.18435 

.03642 

.14176 

-.12007 

-.14361 

.11831 

4.0 

.03741 

.02746 

.07340 

-.08119 

-.07249 

.04035 

5.0 

.01866 

.02634 

.03250 

__-.04345 

r. 03146 

.01299 _ 

6.0 

.01  c8 1 

.01569 

.01345 

-.02058 

-.01299 

.00410 

7.0 

.00672 

.00784 

. .00504 

-.00915 

-.00479 

.00137 

8.0 

.00  28  0 ... 

.00392 

_ .00224  _ 

-.00457 

_. -.0  020  5 . 

..  00068  _ _ 

5.3 

.00112 

.03168 

.00056 

-.00229 

-.00068 

0.00000 

10.0 

.00056 

.00056 

.00056 

-.00114 

0.00000 

0.00000 

Z/Hl 

s 1.01 

X/A  * 4.00 

Y/A 

sx 

SY 

TXY 

SZ  . . 

TXZ  _ 

TYZ  _ . 

1.9 

-.09806 

- .13448 

-.06892 

.07776 

.10805 

.91846 

2.0 

.00672 

-.13000 

-.01905 

.04574 

.02394 

.05881 

3.0 

. 01250 

-.06388 

. .01513  . 

.00457 

. -.01436 

.04377 

4.0 

.02465 

-.01961 

.01737 

-.01029 

-.01710 

.02188 

5.0 

.01431 

- .00280 

•31121 

-.00915 

-.01094 

.00957 

6.0 

. 0 0 72  8 

.00168 

.00560 

-.006  86 

-.00547 

.00342 

7.0 

.00336 

.00168 

.00280 

-.00343 

-.00274 

.00137 

8.0 

.00168 

.00112 

.00112 

-.00229 

-.00137 

.00063 

9.0 

..  .00056 

.00056 _ 

_ .00056 

-.00114 

- .J00JD  68 

JJ..00000 

10.0 

.00056 

.00056 

0.00000 

0.00000 

0.00030 

0 .00039 

Z/Hl 

Y/A 

* 1.01 

sx 

X/A  * 5.00 

SY 

TXY 

SZ 

TXZ 

T YZ 

1.0 

-.03973 

.00734 

-.02017 

-.01144 

.04514 

- .01641  . 

2.0 

-.02353 

- .01961 

-.01961 

.01372 

.02804 

- .00205 

3.0 

-.00560 

-.02129 

-.00784 

.01029 

.01026 

.00615 

4.0 

.00  168 

-,01233 

-.00056 

.00343 

. ..00137 

..  .00615.  . . 

5.0 

.03280 

-.33504 

' .00168 

0.00000 

-.00117 

.00342 

6.0 

. 00168 

- .03168 

.00112 

-.00114 

-.00117 

.00205 

7.0 

.00112 

-.00056 

.00112 

-.00114 

-.00068 

. _ .00068  _ 

8.0 

. 30056 

0 .00030 

•00056 

0.00000 

-.00068 

0 .00000 

9.0 

.03056 

0 .0303C 

0.00000 

0.000  00 

3.00000 

9 .00000 

10.0 

0.03333 

0.00930 

3.00000 

0.00000 

0.00003 

0 .00000 
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TABLE  6 


V 


p = 

.730  NU  = 

.3  00  HI /A 

= .500 

r - * 

E 3/E  1*  3.  000  H3/Hi»  1.00 

E2/E1*  . 

050  H2/H1* 

.02 

* 

1/ HI 

* 1.52  X/A  *_i.00___ 

— — - - 

- 

• - • 

« 

- Y/A 

SX 

SY 

T XY 

sz_ 

TXZ 

tjtz 

1.0 

3.53777 

4 .72429 

5 .18503 

-5.17587 

-8.69858 

.55229 

2.0 

1. 13144 

2.30438 

1.23339 

-2.09955 

-2.36196 

- .32953 

..3.0 

.41530 

.91324 

.32659 

-.80887 

-.64066 

-.,16385. 

4.0 

.15918 

.34913 

.09479 

-.30953 

-.18778 

- .06167 

5.0 

.06  188 

.13307 

.02933 

-.11899 

-.05799 

- .02209 

6.0 

. 02432 

.05115 

.00930 

-.04599 

. -.01933 

-.00736..  .... 

7.0 

. 00966 

.01967 

.00322 

-.017  52 

-.00644 

- .00276 

8.0 

.00358 

.00751 

.00107 

-.00730 

-.00184 

- .00092 

9.0 

.00143 

.00286 

.00036 

-.03292 

-.00092 

0 .00000 

10.0 

.00036 

.00107 

0.C0CG0 

-.00073 

O.COOOO 

0 .00000 

z/Hi  = i.52  y/a  * 2.ao 


Y/A 

sx 

SY 

T XY 

sz 

TXZ 

TTZ 

1.0 

1.63331 

rl. 01053 

.78446 

-.66359. 

_ -1.87227  . 

2.  68781  _ 

2.0 

.75656 

.17313 

•64853 

-.64899 

-1.35311 

.75388 

3.3 

.2  5 46  9 

.24897 

•26542 

-.35187 

-.54032 

.16385 

4.3  „ 

.11447 

.14130 

..09622 

. -.16207. 

..-.19514 

.03314 

5.0 

. 04543 

.06618 

.03398 

-.07008 

-.06904 

.00644 

6.0 

.01824 

.02897 

.01216 

-.02482 

-.02393 

.30092 

7.0 

.00715 

.01216 

.00429 

. . r. 01241 

-.00828  . 

0 .00000  . 

8.0 

.00286 

.00429 

.09143 

-.0  0511 

-.00276 

0 .00000 

9.3 

• 00  107 

.00215 

.00072 

-.00219 

-.00992 

0 .00000 

. ..  1C. 3 

. 03036 

.00072 

.00036 

...  -.00073 

0.00900 

0 .00009 

1 


TABLE  . 6 (Continued) 


Z/Hi 

* 1.52  X/A 

= 3.00 

Y/A 

sx 

sy~ 

TXY 

sz 

TXZ 

T 'll 

1.0 

.14416 

-.42174 

-.01288 

.10147 

-.01657 

.51823 

2.0 

.20211 

-.18565 

.11089 

-.054  75 

26326  _ 

.38752 

3.0 

.11769 

-.02325 

•09050 

-.07665 

-.19330 

.15924 

4.1 

.05580 

.01753 

•04686 

-.05183 

-.09757 

.05431 

5.0 

. 02468 

.01681 

.02075 

..  -.02774.  .. 

.-,04234. 

_ .01749 

ft.  5 

.01073 

.01002 

•00859 

-.01314- 

-.01749 

.00552 

7.0 

.00429 

.00501 

.00322 

-.005  84 

-.00644 

.00184 

S.O 

.00179 

.00250 

_.  00143. 

00292  __ 

_-.00276_ 

„.. 00092  _ 

9.1 

. 00  072 

.00107 

• C0036 

-.00146 

-• C 0092 

C .00000 

10.0 

.00036 

.00036 

.00036 

-.000  73 

0.03000 

0 .00000 

Z/HI 

= 1.52  X/A 

= 4.00 

Y/A 

SX 

SY 

TXY 

SZ 

TXZ  . 

T >Z 

1.1 

-. 06260 

-.08585 

-.04400 

.04964 

•14544 

02485 

2.0 

.03429 

-.08299 

-.01216 

.02920 

.03222 

.07916 

3.1 

.02075 

-.04078 

.00966 

. .00292 

-.01933  . 

.05891 

4.0 

.01574 

-.01252 

.01109 

-.00657 

-.02301 

.0  2946 

5.0 

. 00  894 

-.00179 

.00715 

-.005  84 

-.01473 

.01289 

6.0 

.00165 

.00107. 

.00358 

-.00438 . 

G 0.7.36  __ 

.00460 

7.3 

.00215 

.00107 

.00179 

-.00219 

-.00368 

.00184 

8.0 

.00107 

.00072 

.00072 

-.00146 

-.00184 

.30092 

9.0. 

.00036  _ 

.00036 

_ .00036 

. .-.0  00  73  _ 

..-•0  0092..  _ 

.0  .00000.  ._ 

10.0 

.00036 

.00036 

0.00000 

0.00000 

0.00000 

9 .00000 

Z/Hi  * 1.52  X/A  * 5.00 


- 

- 

-• 

- - 

- -* : ‘ — — 

- — — — -- 

— — — - - — • 

} 

Y/A 

SX 

SY 

TXY 

SZ 

TXZ 

T >Z 

I 

1.9 

-.02540 

.00501 

.-.01288 

-.007  30 

. ...06075 

- .02209 

2.9 

-.01502 

-.01252 

-.01252 

.00876 

.03774 

- .00276 

3.9 

-.00  258 

-.31359 

-.005C1 

.00657 

.01381 

.00828 

4.3 

.00107 

-.00737 

-.00036 

. .00219 

. .00184 

.00828 

I 

5.9 

.09179 

- .00322 

.00107 

0.00000 

-.00184 

.00460 

j 

6.9 

.00107 

-.00197 

.00072 

-.00073 

-.00184 

.00276 

7.0 

. 00072 

-.00036 

_ ..00072 

-.000  73 

-.09092 

.00092 

1 ' 

8.9 

. 00036 

0 .00000 

.00036 

0.000  00 

-.00092 

0 .00000 

% 

9.3 

.09036 

0 .90000 

0.00000 

0.00000 

0.00000 

0 .00000 

1 

1C. 9 

0. 00000 

0 .00009 

0.00000  . 

0.00000 

. 0.00330 

0 .09000 
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TABLE  7 


P * .731  NU  *’.3C0  Hl/A  a 7sOO~ 

E3/E1*  3.000  H3/H1*  1.00  E2/E1*  .050  H2/H1*  .lfl 

Z/Hl  a .50  X/A  * 1.00  _ . _ _ . _ 


Y/A 

. sx 

SY 

TXY 

_ _ sz 

_ TXZ  __  _ 

TYZ  

1.0 

7.10438 

9.48711 

10.41233 

-10.39394 

-3.53990 

.72476 

2.1 

2.27211 

4.62754 

2.47684 

-4.21621 

-.96120 

- .13410 

3.0 

.83  399 

1 .83392 

.65584 

-1.62433 

-.26172. 

T .06668  ... 

4.1 

. 31966 

.70110 

.19036 

-.62158 

-.07642 

- .02510 

5.1 

.12427 

.26722 

.05890 

-.23896 

-.02360 

• .00899 

6.0 

.04885 

.10272  . 

..  .01868 

-.09236 

00787 

_-  .0030  0.  .. 

7.0 

. 01940 

.03951 

•00647 

-.03518 

-.00262 

- .00112 

8.1 

. 01  718 

.01509 

.00216 

-.014  66 

-.00075 

- .00037 

5.0 

.00287 

.01575 

.00072 

-.0  05  86  _ 

. ?.00037_  . 

0 .30000  ... 

10.1 

.00  172 

.00216 

0.00000 

-.00147 

0.00000 

0 .00000 

Z/Hl  » .50  X/A  a 2.00 


Y/A 

sx 

SY 

TXY 

SZ 

TXZ 

TYZ 

1.0 

3.27994 

-2.02931 . 

_.1.57532_ 

-1.33259  . 

__r.76192  ... 

.1  .09381 

2.0 

1.51929 

.34768 

1.30235 

-1.30327 

-.55065 

.30679 

3.0 

.51146 

.49996 

.53301 

-.70661 

-.21989 

.06668 

4.1 

.22987 

•28374 

. .19323 

_ - . 32545 

-.07941  .. 

•01349_. 

5.1 

.09123 

.13289 

.06824 

-.140 74 

-.02899 

.00262 

6.0 

•03664 

.05819 

.02442 

-.04984 

-.00974 

.00037 

7.0 

.01437 

. .02442 

.00862. 

r.  024  92  .. 

..-•00337. ... 

0 .00000  ... 

8.0 

.03575 

•00862 

.00287 

-.010  26 

-.00112 

0 .00000 

9.0 

.00  216 

.00431 

.00144 

-.00440 

-.00337 

0 .00000 

10.3 

.00072 

.00144  . 

•00072 

. -..00147 

0.00000 

0 .00000. 

TABLE  7 (Continued) 


Z/Hl  = .50  X/A  = 3.00 


Y/A 

sx 

SY 

TXY 

sz 

TXZ 

T YZ 

1.9 

.28949 

>.84692 

-.02586 

.20377 

-.00674 

21090 

2.0 

• 40586 

. - .37282  . 

.22269 

..._-.10995  _ 

_-.107_13__ 

_ .15770.  _ 

3.0 

.23633 

- . 04669 

.18174 

-.15393 

-.07866 

.0  6480 

4.0 

.11206 

.03520 

.09410 

-.104  09 

-.03971 

.02210 

5.0 

.04  95  7 

.03376 

.04166 

-.05571 

-.01723  _ 

007,12 

6.0 

.02155 

.02011 

.01724 

-.026  39 

-.00712 

.00225 

7.0 

•00862 

•01006 

.00647 

-.01173 

-.C0262 

.00075 

8.0 

.00359 

.00503 

. 00  287 

.-..00586  . 

00112  . 

.000  37 

9.0 

.00144 

.00216 

.00072 

-.00293 

-.00037 

0 .00030 

10.0 

.00072 

.00072 

.00072 

-.00147 

0.00300 

0 .00000 

Z/Hl 

* .50  : 

X/A  = 4.00 

Y/A 

sx 

SY 

TXY 

sz 

TXZ 

* T.YZ 

1.0 

-.12571 

-.17243 

-.08836 

.09969 

.05919 

.01011 

2.0 

. 00862 

-.16665 

-.02442 

•05864 

•01311 

.0  3221 

3.0 

. 04  166 

- .08189 

.01940 

....  .00586 

r. 00737 

.02397  . 

4.0 

.03161 

-.02514 

.02227 

-.01319 

-.00936 

.01199 

5.0 

.01796 

-.00359 

.01437 

-.01173 

-.00599 

.00524 

6.0 

_ .00934 

.30216  . 

.00718 

.-..0  08  80 

-.00300, 

.00187 

7.3 

. 00  431 

.00216 

.00359 

-.0  0440 

-.00150 

.00075 

8.0 

.00216 

.00144 

.00144 

-.00293 

-.00075 

.00037 

9.0 

. 00  072 

...00072  _ 

00072 

r.  00147.. 

.-.00337. 

__  0 .00000  . 

10.0 

. 00072 

.00072 

0.00000 

0.0  00  00 

0.00000 

0 .00000 

Z/Hl 

* .50 

X/A  = 5.00 

Y/A 

SX 

SY 

TXY 

sz 

TXZ 

T >Z 

1.0 

-.05103 

.01006 

-.02586 

-.01466 

. Q 2472 

-..00899  ... 

2.0 

-.03017 

-.02514 

-.02514 

. 0 17  59 

.01536 

- .00112 

3.0 

-.00718 

-.02730 

-.01006 

.01319 

.00562 

.00337 

4.3 

.00216 

-.01580 

-.00072 

• 00440 

•00075 

.00337 

5.0 

.00  359 

-.00647 

.00216 

0.00000 

-.00075 

30187 

f.O 

.00216 

-.00216 

.00144 

-.00147 

-.00075 

.00112 

7.0 

.00144 

-.00072 

.00144 

-.00147 

-.000  37 

.00037 

8.  0 

.00072 

C. 00003 

.30072 

0. 0 00  CO 

-.00037 

0 .00000 

9.0 

. 03072 

0.00030 

0.00000 

0.00000 

0.00000 

0 .00000 

10.0 

0.00000 

0 .03030 

0.00000 

0.00000 

0.00000 

0 .00000  _ 

67 


TABLE  8 


p = 

.700  NU  * 

.3  00  Hl/A 

* .500 

E3/E 

1*  3.000  H3/H 1*  1.00 

E2/E1*  , 

050  H2/H1* 

.10 

Z/Ml 

*1.05  . X/A  a 1.30 

• 

y/a 

_S  X . . . 

_ SY  _ . .. 

T.Xr.  _ 

sz 

. TXZ* 

. T YZ__  . _ 

1.0 

5.61039 

7.49205 

0.22271 

-0.20010 

-6.92279 

.43954 

2.0 

1.79430 

3.65441 

1.95590 

-3.32950 

-1.07977 

- .26226 

3.0 

.65061 

.1.44026 

.51793 

-1.20275 

_ r. 50907 

.13040 

4.9 

.25  244 

.55366 

.15033 

-.490  07 

-.14944 

- .04900 

5.0 

. 09014 

.21103 

.04652 

-.10071 

-.04615 

- .01750 

6.0 

.0305  7 

.00112 

.01475 

-.07294 

-.01530 

.-.005  06 

7i0 

.01512 

.03120 

' .00511 

-.02779 

-.00513 

- .03220 

6.0 

.00567 

.01191 

.00170 

-.01150 

-.00147 

- .00073 

9.0  . 

.00227 

.00454 

.00057 

.-,00463 

-.00073 

0 .00000  _ 

10.3 

.00057 

.00170 

0.00000 

-.00116 

0.00000 

0 .00000 

Z/Hl 

* 1.05  X/A 

* 2.00 

Y/A 

SX 

SY 

TXY 

SZ 

TXZ 

T>Z 

1.0 

2.59020  -1 

.60256 

1.24404 

-1.05236 

_?1.49.0  05 

-2.13910 

2.3 

1.19900 

.27456 

1.02040 

-1.02921 

-1.07606 

.59997 

3.0 

.40390 

.39403 

.42092 

-.55002 

-.43002 

.13040 

4.1 

. 16  153 

.22406 

.15260 

-“.25701 

. -.15530 

. .02637_ 

5.0 

.37204 

.10495 

.05309 

-.11114 

-.05494 

.00513 

6.0 

.02093 

.04595 

.01929 

-.03936 

-.01905 

.00073 

7.0 

.01  135 

.01929 

-.00601 

019  60. 

..  _r. 00659 

0 .00000  . 

0.0 

.00454 

.00601 

.00227 

-.00010 

-.00220 

0 .00000 

9.0 

• 00170 

.00340 

.00113 

-.00347 

-.09073 

0 .00000 

10.0 

.00057 

.00113 

. .00057 

-j.  0 01 16 

_ 0.00000 

_C  .00000  . . 

I 


TABLE  8 (Continued) 


Z/Hl 

* 1.05  X/A 

= 3.00 

Y/A 

sx 

SY 

TXY 

sz 

TXZ 

T *Z 

1.0 

.22861 

•66682 

-.02042 

.16092 

-.01319 

■41244 

2.9 

.32051 

•29442 

...  .17586 

-.066SJ 

r . 20952.  . 

_ .30841 

3.0 

.18663 

•03687 

.14352 

-.12156 

-.15384 

.12673 

A. 9 

.08650 

.02760 

•07431 

-.08220 

-.07765 

.04322 

5.0 

.03914 

.02666 

...  ..03290 

-.04399. 

0 3370 

01392 

6.0 

.01702 

.01568 

.01361 

-.02084 

-•C1392 

.00440 

7.9 

.00661 

.00794 

•00511 

-.00926 

-.00513 

.00147 

6.0 

.00284 

.00397 

_ .00227 

-.0  0463 

0 0 220.  . 

.00073 

5.0 

.00  113 

.00170 

.00057 

-.00232 

-.00073 

0 .00000 

10.0 

.00057 

.00057 

•00057 

-.00116 

0.00000 

0 .00000 

Z/Hl 

= 1.05 

X/A  * 4.00 

Y/A 

sx 

SY  _ 

. TXY  . _ 

sz 

TXZ  . . . 

T_>Z 

1.9 

-.09927 

- .13615 

-.06978 

.07872 

.11575 

.01978 

2.0 

.00  P81 

-.13161 

-.01929 

• 04631 

.02564 

.06300 

3.9 

. 03290 

-.06467 

.01532 

.00463 

-.01538. 

. .0468.9 

4.9 

.02496 

-.01985 

.01759 

- . 0 10  42 

-.01831 

.02344 

5.3 

.01418 

- .00284 

.01135 

-.00926 

-.01172 

.01026 

6.0 

. 00737 

.00170 

.00567 

-.0  0695 

-.00586 

.00366 

7.3 

.00340 

.00170 

.00294 

-.00347 

-.00293 

.00147 

8.9 

• 00  170 

.00113 

.00113 

-.00232 

-.00147 

.00073 

9..0 

..00057 

_ . .00057 

. .00057 

_?•  0 0116 

-.00073 

0 .00000 

10.0 

.00057 

.00057 

0.00000 

0.000  00 

0.00000 

0 .00000 

Z/Hl 

* 1.05 

X/A  * 5.00 

V/A 

SX 

SY 

TXY 

SZ 

TXZ 

T >Z 

1.0 

-.04028 

.00794 

-.02042 

-.01158 

•C4835 

- .01758 

2.0 

-.02383 

- .01985 

-.01985 

.01389 

.03004 

-.00220 

3.3 

-.00567 

-.02156 

-.0  0794 

.01042 

.01099 

• 006S9 

4.9 

.00170 

-.01248 

-.00057 

.0  0347. 

.00147  . 

.00659 

5.9 

.00284 

-.00511 

.00170 

C. 000  00 

-.00147 

.00366 

6.0 

.00170 

-.00170 

.00113 

-.00116 

-.00147 

.00220 

7.0 

.00  113 

-.00057 

_ .00113.  ... 

T.  00116 

-.00073.. 

_ .00073 

8.0 

.00057 

0.00300 

•00057 

0.00000 

-.00073 

0 .00000 

9.9 

.00057 

0 .00000  ■ 

0.00000 

0.00000 

0.00000 

0 .00000 

10.3 

0.00000 

0.00000 

0.00000  . . 

0.00000 

0.00000. 

0 .00000 

TABLE  9 


P * . 

700  NU  » . 

300  HI /A 

* .500 

r 

t 

E3/E1 

* 3.000  H3/H 1*  1.00 

E2/E1*  . 

050  H2/H1* 

.10 

1 i 

. 7/Hl 

* 1.60  X/A 

» 1.00 

- 

— - - 

- - — 

V/A  _ 

_ sx 

SV  . 

..TXY 

sz 

_ T.X7  . _ 

T >Z  . 

1.3 

3.29502 

4.40014 

4.82926 

-4.82072 

-9.29216 

.58998 

2.3 

1.05381 

2.14626 

1.14876 

-1.95549 

-2.52314 

- .35202 

V 

3.0 

.36661 

.85058  .. 

. .30418 

7.53  37 

. -.68437 

- .17503 

4.0 

.14826 

.32517 

.08829 

-.28829 

-.20059 

- .06588 

5.0 

.05  764 

.12394 

.02732 

-.110  83 

-.06195 

- .02360 

6.3  . . 

.02266  . 

.04764  .. 

..  .00  866 

- .042-84 

_-. 02065 

-..00787. 

7.0 

• 33900 

.01832 

.00300 

-.01632 

r. 00688 

- .00295 

P.3 

.00333 

.30700 

.00100 

-.00680 

-.00197 

- .00098 

_9  • 0 

• 00  13  3.  _ 

.00267 

.00033  . 

..-.00272 

-.00398 

0 .00000  . 

* J 

10.0 

.00033 

.00100 

0.00000 

-.00068 

0.00000 

0 .00000 

Z/Hl 

* 1.60  X/A 

* 2.00 

1 

Y/A 

SX 

SY 

TXY 

SZ 

TXZ 

T YZ 

1 - . 

.1.3 

1.52124 

-.94120  _ 

. .73064  . 

-.61606  _ 

,r2. 00  03  3 

2..87123 

2.0 

.70465 

.16125 

•60403 

-.60446 

-1.44545 

.80532 

3.3 

.23721 

.23186 

.24721 

-.32773 

-.57720 

.17503 

4.0 

.10661  . 

.13160  .. 

.08962, 

-.15095  __ 

-.20846  _ 

..0  3540  

T~~  r *• 

5.0 

.04  231 

.06164 

.03165 

-.065  27 

-.07375 

.00688 

6.0 

.01699 

.02699 

.01133 

-.02312 

-.02557 

.00098 

7.3 

.00666  . 

•01133  . 

-.00400 

. - . 01156  . . 

r. 00885 

0.00000  

6.0 

.00267 

.00400 

.00133 

-.0  04  76 

-.0  0295 

0.00000 

9.0 

• 00  100 

.03200 

.00067 

-.00204 

-.00096 

0.00000 

i - • 

10.0 

•00033 

.00067 

.00033 

. -.00068  . 

,_0. 00000  . 

0.00000 

TABLE  9 (Continued) 


Z/Hl 

*1.60  X/A 

* 3.00 

Y/ A 

sx 

sf  “ 

TXY 

SZ 

TXZ 

TYZ 

1.0 

.13427 

•39280 

-.01199 

.09451 

-.01770 

.55360 

2.0 

.18  624 

.17291 

_ .10  328.. 

- -.050 99 _ 

-.28122 

..41397, 

3.3 

.10961 

.02166 

.08429 

-.07139 

-.20649 

.17011 

4.0 

. 05197 

.01633 

.04364 

-.04828 

-.10423 

.05801 

5.0 

_ _. 02299 

.01566 

.01932 

_-.025  84_ 

=•04523. 

,01868 

6.0 

.01003 

.00933 

.00800 

-.01224 

-.01868 

.00590 

7.0 

.03403 

•00466 

•00300 

-.00544 

-.00688 

.00197 

8.0 

.03  167  _ . 

.00233 

.00133 

. .-.00272 

-.00295 

.00098 

9.0 

.00067 

.00100 

•00033 

-.00136 

-.00098 

0 .00000 

10.0 

. 00  033 

.000  33 

.00033 

-.00068 

0.00000 

0 .00000 

Z/Hl 

* 1.60 

X/A  * 4.00 

Y/A 

sx 

. SY  . ... 

TXY.  

SZ  _ 

TXZ. 

TYZ.  . 

1.1 

-.05830 

-.07996 

-.04098 

.04624 

.15536 

02655 

2.0 

.00  <<00 

-.07729 

-.01133 

•02720 

.03442 

.08456 

3.1 

. .01932 

-.03798 

•QC9C0  . 

.00272. 

r. 02065  . 

...  .06293 

4.0 

.01  t66 

-.01166 

•01033 

-.00612 

-.02458 

.03147 

5.0 

. 00833 

-.00167 

.00666 

-.00544 

-.01573 

.01377 

6.0 

. 00  t3  3 

.00100 

.00333 

-.0  04  08 

r. 00787 

. .00492 

7.0 

. 00200 

.00100 

.00167 

-.002  04 

-.CC393 

.00197 

8.0 

.00  100 

.00067 

.00067 

-.00136 

-.00197 

.00098 

.9.1  .. 

_ .03033 

.80033  _ 

• 00033  . , 

-.00068 .... 

r.00098. 

JL.00000, 

10.1 

.00  033 

•00013 

0.00000 

0.000  00 

0.00000 

0 .00000 

Z/Hl 

* 1.60  X/A  * 5.00 

Y/A 

SX 

sv 

TXY 

SZ 

TXZ 

T YZ 

1.0 

-.02365 

.00466 

-.01199 

. -.00680 

• Q6490 

- 02360 

2.0 

-.01399 

-.01166 

-.01166 

.00816 

.04032 

- .00295 

3.9 

-.00333 

-.01266 

-.00466 

.00612 

•C1475 

.00885 

4.0 

.,00100 

-.00733 

-.00933 

,.Q0204_ 

„ .00197 

_. 00885 

5.9 

.00167 

-.00300 

.00100 

0.000  00 

-.00197 

.00492 

6.0 

.00100 

- . 00100 

.00067 

-.00068 

-.00197 

.00295 

7.1 

.09  067 

-.00033 

,00067 

-.00068. 

. r.00098 

.00098 

8.0 

.00033 

0.00000 

.00033 

0.000  00 

-.09098 

0 .00000 

9.9 

.00033 

0.00000 

0.00000 

0.00000 

0.00000 

0 .00000 

10.0 

0.09000 

0 .00009 

0.00000 

...  0. 00000 

0.00000 

0 .00000 
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DISTRIBUTION  LIST 
Govarnnant  Activities 

No.  of 
Copies 


NAVAIRSTSCGM,  AIR-954  (2  for  retention),  2 for  AIR-530, 

1 for  AIR-320B,  AIR-520320,  AIR-5302,  AIR-53021,  AIR-530215).  9 

APFDL,  WPAFB,  OH  45433 

(Attn  : FB/Mr.  P.  A.  Panaley) 2 

(Attn:  FBC/Hr.  C.  Wallace) 1 

(Attn:  PBC/Mr.  E.  E.  Zink)  .................  1 

AFML,  WPA71,  OR  45433  ^ 

(Attn:  LAM  (Technical  Library)) 1 

(Attn:  LT-l/Mr.  W.  R.  Johnston) 1 

(Attn:  LIF/Kr.  T.  Cordell) 1 

(Attn:  FBSC/Mr.  L.  Kelly) 1 

(Attn:  MAC /Mr.  G.  P.  Peterson) 1 

(Attn:  MXlA/Mr.  F.  J.  Fechefc) 1 

(Attn:  MBC/Kr.  T.  6.  Ralnhard,  Jr.) % 1 

AFOSR,  Washington,  D.C.  20333 

(Attn:  Mr.  J.  Ponte  rants) 1 

DOC 12 

FAA,  Alrfraaas  Branch,  PS- 120,  Washington,  O.C.  20553 
. (Attn:  Mr.  J.  Dougherty)  ••••••••••..•••••.  1 

NASA  (ADM),  Washington,  D.C.  20546 

(Attn:  Seeratary).  . 1 

RASA,  George  C.  Marshall  Space  Fllsht  Center,  Huntsville, 

AL  35812 

(Attn:  Stf-ASTN-ES/Mr.  I.  I.  Rnglar)  ‘l 

(Attn:  SAB -ASTN-M/Mr  j R.  Schirlnghansr) 1 

(Attn:  8tf-ASlM-ltlM/Dr.  J.  M.  Stuckey)  ...........  1 

NASA,  Langley  Research  Center,  Haavton,  VA  23365 

(Attn:  Mr.  J.  P.  Peterson,  Mr.  R.  Pride,  and  Dr.  M.  Card).  3 
■ASA,  Lewis  Research  Center,  Cleveland,  OH  44153 

(Attn:  Technical  Library,  and  M.  Hershberg).  . 2 

MAVPGSCHL,  Monterey,  CA  95940 

(Attn:  Prof.  R.  Ball,  Prof.  M.  H.  Bask) 2 

KAVS1ASY8C0M,  Washington,  D.C.  20362 

(Attn:  Code  035,  Mr.C.  Pohler).  •••••••  1 

KAVSEC,  Ryattsvllle,  MD  20782 

(Attn:  Code  6101E03,  Mr.  W.  Craner).  ......  1 

MAVMZPRAMDC8M,  Bethesda,  » 20034 

(Attn:  Code  173.2,  Mr.  W.  P.  Couch) 1 

NAVSHIPIAXDCSN,  Annapolis,  MD  21402 

(Attn:  Code  2870,  Mr.  H.  Edelstain).  ••••••••••••  1 

ML,  White  Oak,  MD  20910 

(Attn:  Mr.  F.  R.  Barnet) 1 

ML,  Washington,  D.C.  20375 

(Attn:  Dr.  X.  Nolock) . 1 

OMR,  Washington,  D.C.  203 
(Attn:  Dr.  I.  Pa  none) 


i 

i 

i 

1 
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Csranant  Activities  (Cont.) 

PLASTIC,  Picatinny  Arsenal,  Dover*  NJ  07801 

(Attn:  Librarian,  Bldg.  176,  SARPA-FR-M-D  and  Mr.  H.  Pebly).  • 
Scientific  6 Technical  Info  nation  facility.  Collage  Pack,  MD 

(Attn:  VASA  lepras entatlva) . . . . 

USAAVMATLAB,  Pert  Eustis,  PA  23603 

(Attn:  Mr.  1.  Bereeford)  •......•  . 

OSAMATRESAC,  Vatertovn,  HA 

(Attn:  Dr.  E.  Lenoe)  

USAKESOFC,  Durham,  VC  27701  


2 

1 

1 

% 

1 

1 


i 


Hon-Cov.nnrae  A$racl«s  . 


Avco  Aero  Structures  Division,  Nashville,  IN  37202  . 

(Attn:  Mr.  W.  Ottenville).  1 

Battalia  Columbus  Laboratories,  Metals  and  Ceramics  Inforaatlon 

Center,  505  Xing  Avenue,  OH  43201.  .•••••••••••••  1 

Bell  Aerospace  Coafany,  Buffalo,  NT  14240 

(Attn:  Zone  1*85,  Mr.  P.  M.  Anthony)  • •.•••%••••••  1 

Bell  Helicopter  Company,  Port  Worth,  TX  76100 

(Attn:  Mr.  Charles  Harvey)  • •••••.•••••.  ......  1 

Bend lx  Products  Aerospace  Division,  South  Bend,  IN  46619 

(Attn:  Mr.  R.  V.  Cervelll)  • . • ••••••«•  1 

Boeing  Aerospace  Company,  P.0.  Box  3999,  Seattle,  KA  98124 

(Attn:  Code  206,  Mr.  R.  E.  Horton)  .........  1 

J Boeing  Company,  Renton,  Washington  98055 

(Attn:  Dr.  X.  June).  ••••••••••••.•••»••••  1 

Boeing  Company,  Vertol  Division,  Phils.,  PA  19142  * 

(Attn:  Mr.  1.  L.  Pinckney,  Mr.  D.  Hoffstedt)  •••••••••*2 

Boeing  Company,  Wichita,  XS  67210 

(Attn:  Mr.  Reneau/MS  16-39)  • ••••» • . • 1 

Cabot  Corporation,  Billerica  Research  Center,  Billerica,  MA 

01821 1 

Draxel  University,  Phils.,  PA  19104 

(Attn:  Dr.  P.  C.  Chou)  ••*...•  ••••••••  1 

E.I.  DuPont  Company,  Wilmington,  DE  19898 

(Attn:  Dr.  Carl  Zveben)  Bldg.  262/Room  316  .•••••••••  1 

falrehlld  Industries,  Hagerstown,  MD  21740 

(Attn:  Mr.  D.  Ruck).  1 

Georgia  Institute  of  Technology,  Atlanta,  GA 

(Attn:  Prof.  W.  H.  Horton)  1 

General  Dynamics /Convair,  San  Diego,  CA  92138 

(Attn:  Mr.  D.  R.  Dunbar,  W.  G.  Schcck)  2 

General  Dynamics,  Port  Worth,  TX  76101 

(Attn:  Mr.  P.  D.  Shoe  key,  Dept.  23,  Mall  Zone  P-46). 1 

General  Electric  Company,  Pbila.,  PA  19101 

(Attn:  Mr.  L.  McCrelght)  ...•••••••.....••.»  1 

. .Great  Lakes  Carbon  Corp.,  N.T.,  NT  10017 

(Attn:  Mr.  W.  R.  Bonn,  )fcr«,  Marftay  Development)  .......  1 

C ramson  Aerospace  Corporation,  Bethpago, 

I (Attn:  Mr.  1.  Hadcoek,  Mr.  S.  Destin) 
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Mon-Gov* rumen  t Agencies  (Cone.) 

Hercules  Powder  Company,  Inc. , Cumberland , HD  21501 

(Atta:  Mr.  D.  Hug) 1 

I.  I.  Thompson  Fiber  Glass  Company,  Cardans,  CA  90249 

(Atta:  Mr.  M.  Jfyers) 1 

XXX  lose arch  Institute,  Chicago,  1L  60616 

(Atta:  Mr.  K.  Hofar)  •••••••••••••••••••••  1 

J.  P.  Stevens  A Co.,  Inc.,  H.Y.,  NY  10036 

(Attn:  Mr.  H.  I.  Shulocfc) 1 

Fawn  Aircraft  Corporation,  Bloomfield,  CT  06002 

(Atta:  TSch.  Library) 1 

Lehigh  University,  Bethlehem,  PA  18015 

(Attn:  Dr.  G.  C.  Slh) 1 

Lockheed -California  Company,  Burbank,  CA  91520 

(Attn:  Mr.  B.  K.  Walker,  R.  L.  Vaughn)  ..J1’ ••••••.••  2 

Lockheed -Georgia  Company,  Marietta,  GA 

(Atta:  Advanced  Composites  Information  Canter,  Dept.  72-14, 

Zone  42).  •••••••  •••••••••••  1 

LTV  Aerospace  Corporation,  Dallas,  T X 75222 

(Atta:  Mr.  0.  E.  Dhonau/2 -53442 , C.  R.  Foreman) 2 

Martin  Company,  Baltimore,  MD  21203 

(Attn:  Mr.  J.  E.  Favicon)  1 

Materials  Sciences  Corp. , Blue  Bell,  PA  19422  ...«••••••  1 

McDonnell  Douglas  Corporation,  St.  Louis,  MO  63166 

(Attn:  Mr.  B.  C.  Goran,  0.  B.  Me Bee,  C.  Stenberg).  ......  3 

McDonnell  Douglas  Corporation,  Long  Beach,  CA  90801 

(Attn:  H.  C.  Schjulderup,  G.  Lehman)  2 

Minnesota  Mining  and  Manufacturing  Company,  St.  Paul,  M 55104 

(Attn:  Mr.  W.  Davis) 1 

Northrop  Aircraft  Corp.,  Moral r Dlv.,  Hawthorne,  CA  90250 

(Attn:  Mr.  B.  D.  Hayes,  J.  V.  Noyes,  R.  C.  Isemann).  . • • . . 3 

Rockwell  International,  Columbus,  OH  43216 

(Attn:  Mr.  0.  G.  Acker,  K.  Clayton) 2 

Rockwell  International,  Los  Angeles,  CA  90053 

(Attn:  Dr.  L.  Lackman)  ....................  1 

Rockwell  International,  Telsa,  OK  74151 

(Attn:  Mr.  E.  Sanders,  Mr.  J.  R.  Powell)  •••••••••..  2 

Omens  Corning  Fiberglass,  Granville,  OH  43023 

(Attn:  Mr.  D.  Mattes) 1 

Rohr  Corporation,  Riverside,  CA  92503 

(Attn:  Dr.  F.  Biel  and  Mr.  R.  Elkin) 2 

Ryan  Aeronautical  Cosy  any,  San  Diego,  CA  92112 

(Attn:  Mr.  R.  Long) 1 

Sikorsky  Aircraft,  Stratford,  CT  06497 

(Atta:  Mr.  J.  Ray) • l 

University  of  Oklahoma,  Norman,  OK  93069 

(Atta:  Dr.  G.  M.  Nordby)  •••••••••••••••••••  1 

Union  Carbide  Corporation,  Cleveland,  OH  44101 

(Attn:  Dr.  H,  F.  Volk) 1 


